Photoassociative creation of ultracold heteronuclear ^Li^°K* molecules 
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We investigate the formation of weakly bound, electronically excited, heteronuclear ''Li*°K* 
molecules by single-photon photoassociation in a magneto-optical trap. We performed trap loss 
spectroscopy within a range of 325 GI^z below the Li(25'i/2)+K(4P3/2) and Li(25'i/2)+K(4Pi/2) 
asymptotic states and observed more than 60 resonances, which we identify as rovibrational lev- 
els of 7 of 8 attractive long-range molecular potentials. The long-range dispersion coefficients and 
rotational constants are derived. We find large molecule formation rates of up to ~ 3.5 x 10^s~^, 
which are shown to be comparable to those for homonuclear ^"K^. Using a theoretical model we 
infer decay rates to the deeply bound electronic ground-state vibrational level X^E"'"(w' = 3) of 
~ 5 X 10*s~^. Our results pave the way for the production of ultracold bosonic ground-state ®Li*°K 
molecules which exhibit a large intrinsic permanent electric dipole moment. 



I. INTRODUCTION 



The recent realization of gases of ultracold polar 
molecules in their rovibrational ground state [H Q has 

fened a new frontier in atomic and molecular physics 0, 
Due to their long-range anisotropic dipole-dipole in- 
teractions and the possibility to trap and manipulate 
these molecules with external electric fields, they offer 
fascinating prospects for the realization of new forms of 



quantum matter 
mation processing 



6^, Applications to quantum infor- 
precision measurements 
and ultracold chemistry [3| have been proposed. 

The heteronuclear alkali dimer LiK is an excellent can- 
didate for these studies. It has a large dipole moment 
of 3.6 D 11] in its singlet rovibrational ground state 
and both of its constituents, Li and K, possess stable 
fermionic and bosonic isotopes with which dipolar gases 
of different quantum statistics can be realized. 

While atoms are routinely laser cooled to ultra- 
cold temperatures, the complex internal structure of 
molecules makes this direct method difficult (although 
possible fl^l). So far the most efficient way to produce ul- 
tracold molecules has been to associate pre-cooled atoms. 
Two techniques have been established, namely magnet- 
ically tunable Feshbach resonances and photoassocia- 
tion. Feshbach resonances allow the production of vibra- 
tionally excited molecules in the electronic ground state. 
In this way, ultracold heteronuclear ^Li'^^K molecules 
could recently be produced [H, [l4| . A combination of 
Feshbach resonances with a multi-photon state transfer 
may give access to the collisionaUy stable rovibrational 
ground state [l|, [l^ . Photoassociation can directly give 
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access to this state either via single-photon photoassoci- 
ation and subsequent s pon taneous decay or by multi- 
color photoassociation |16| . 

In this letter we report on the production of ultracold 
heteronuclear excited ^Li^^K* molecules by single-photon 
photoassociation (PA) in a dual-species magneto-optical 
trap (MOT). We detect the molecule creation by a loss 
in the number of trapped atoms, which results from the 
molecules' spontaneous decay into either a pair of free 
untrapped atoms or a bound ground-state molecule. 

Heteronuclear PA has so far been demonstrated for 
RbCs* [l3, KRb* [11, NaCs* [H, LiCs* and 
YbRb* [201 . As compared to homonuclear molecules, the 
PA rate for heteronuclear molecules is typically smaller 
due to the different range of the excited-state potentials. 
Whereas two identical atoms in their first excited state 
interact via the resonant dipole interaction at long range 
(with potential V{R) oc —Cs/R^), two atoms of differ- 
ent species interact via the van der Waals interaction 
{V(R) oc —Ce/R^), leading, for the heteronuclear case, 
to molecule formation at much shorter distances where 
fewer atom pairs are available. Besides, it has been ar- 
gued that among the heteronuclear dimers, LiK* would 
be particularly difficult to photoassociate due to its small 
reduced mass and Cg coefficients, which should lead to 
small PA rates of e.g. two orders of magnitude less than 
for the heavier dimers RbCs* and KRb* [21]. However, 
the PA rates we observe in our experiment are similar 
to those of the comparable experiment with RbCs* [13] 
and those found for homonuclear Kj. Our theoretical 
calculations are able to predict the large rates observed. 

We perform PA spectroscopy in order to determine the 
long-range part of the excited-state molecular potentials. 
Previously, several molecular potentials of LiK have been 
determined by molecular |22h24| and Feshbach resonance 
spectroscopy fl4, 2^. Our measurements give access to 
previously undetermined spectroscopic data of high pre- 
cision and will allow the derivation of more precise molec- 
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ular potential curves facilitating the search for efficient 
pathways to produce LiK molecules in the rovibrational 
ground state. 

Figure [Ija) shows the molecular potentials dissociat- 
ing to the three lowest electronic asymptotes 25+45', 
25+4F and 2P+45 of the LiK molecule. They have 
been calculated as described in ref. a nd connected 
to the asymptotic form given in ref. [2a] at large dis- 
tances {R > 40 oq). Note that as usual in alkali dimers, 
a strong spin-orbit coupling is expected between the l^II 
and the 2^E+ states due to the crossing of their potential 
curves around 7.5 oq. Relevant for our experiment are the 
eight Hund's case c potential curves dissociating to the 
25i/2+4Pi/2,3/2 asymptotes. Figure [TJb) displays their 
long-range part, which is obtained by diagonalizing the 
atomic spin-orbit operator in the subspace restricted to 
the Hund's case a states correlated to 25-I-4P, for each 
of the symmetries fi*^ = 0"'",0~,1,2 (where ft denotes 
the quantum number of the projection of the total elec- 
tronic angular momentum on the molecular axis and a 
the parity of the electronic wave function through a sym- 
metry with respect to a plane containing the molecular 
axis). These potentials are all attractive at long range, 
whereas the curves which dissociate to the asymptotes 
2Pi/2,3/2+45i/2 are all repulsive [2l|. For the relevant 
asymptotes the dispersion coefficients Cg assume only 
three different values due to the small atomic fine struc- 
ture of the Li atom [27| . They have been calculated theo- 
retically [2^, [131 they are determined experimentally 
in this work. 

II. EXPERIMENTAL SETUP 

The ^Li'^^K* molecules are created by a PA beam 
which is superimposed with the atoms trapped in the 
dual-species MOT. The MOT is continuously loaded 
from a Zeeman slower for ^Li and a 2D-M0T for ^°K, as 
described in ref. ^2^ . We record PA spectra by scanning 
the PA beam in frequency, simultaneously recording the 
steady-state atom number of each species via the emitted 
trap fluorescence. The signature of ^Li'^^K* formation is 
a decrease of both the ^Li and the fluorescence. The 
PA laser is scanned red detuned with respect to one of 
the atomic transitions of '*°K (see fig. [Ijb)) and has no 
effect on a single-species ^Li-MOT. The ^Li fluorescence 
signal thus represents a pure heteronuclear PA spectrum, 
whereas the ^'^K fluorescence signal represents the sum 
of a heteronuclear (^Li^°K*) and homonuclear ('"'Kj) PA 
spectrum. The frequency of the PA laser is recorded by a 
wavelength meter (High Finesse, ref. WS-6) with an abso- 
lute accuracy of ±250 MHz. Additionally, a Fabry- Perot 
interferometer is used to verify the laser's single-mode 
operation. 

The PA light is derived from a homemade diode laser- 
tapered amplifier system. It has a wavelength of 767 nm 
and a power of 660 mW at the output of a single- 
mode polarization-maintaining fiber. It is collimated and 
passes four times through the center of the MOT with a 
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FIG. 1: (a) Molecular potentials of the LiK molecule for short 
interatomic separations R (ao = 0.0529177 nm). The upward 
arrow represents the energy delivered by the PA laser and 
the downward the spontaneous decay to electronic ground- 
state molecules. The vibrational state X^E"'"(v' = 3) shown 
in the figure has a favorable overlap with the addressed ex- 
cited states due to spatially coincident classical inner turning 
points, (b) Detailed view of the excited-state potentials, la- 
beled by their Hund's case c quantum numbers and an 
additional classification ("up/down") for unambiguous dis- 
tinction. At short range, each of these potentials approaches 
one of those shown in fig. (a) as illustrated in ref. [281 ]. 



total peak intensity of ~ lOOW/cm^. The beam diam- 
eter of 2.2 mm (1/e^) was chosen to match the size of 
the ^Li-MOT. Using the feed-forward technique [13], the 
laser's mode hop free continuous tuning range extends 
over ~ 35 GHz. 

For optimum experimental conditions, the PA-induced 
trap loss needs to be maximized and all other intrinsic 
losses that compete with it minimized 28]. Besides, the 
frequency of the PA beam needs to be scanned slowly 
enough 15MHz/s) to allow the trap loss to reach 
a quasi-steady state. To achieve these conditions the 
^Li-MOT is reduced to a small atom number and vol- 
ume (by lowering the loading rate) and placed at the 
center of the larger '"^K-MOT. Further, light-induced 
cold collisions are reduced by using small intensities for 
the MOT cooling and repuniping light {1^^^^ ~ 1.5/,'"' 

orK 

spectively). The detunings are Ai/^^^^ 



Li 

cool ' ^ -^■'-'-'sat' 

104^t, ^rep ~ S/.'^t per beam, re- 
"01 = A^.^^ ~ -3r, 



-4r and the axial magnetic field gradi- 
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ent is 20 G/cm. These parameters result in A^Li '-^ 5 x 10* 
and A'k ~ 2.5 x 10^ trapped atoms with central atomic 
densities of nu ^ 7 x 10^" cm~'^ and tik ~ 5 x 10^° cm^'^ 
and temperatures of ^ 1.2 mK and Tk ~ 300 fiK, 
respectively. At these temperatures only heteronuclear 
collisions of s- and p-wave character (i.e., £ — 0,1, where 
£ is the rotational angular momentum of the atom pair) 
reach sufficiently short internuclear distances to allow for 
PA (the height of the d-wave rotational barrier being 
13.4 mK). If J is the total angular momentum quantum 
number of the atom pair (including electronic angular 
momentum and rotation), molecule formation is thus re- 
stricted to the rotational levels J = 0, 1,2 for electronic 
states with f2 = 0, J = 1, 2, 3 for £7 = 1 and J = 2, 3 for 

n = 2. 

III. RESULTS 

Figure [2] (a-d) shows a compilation of our recorded 
spectroscopic data. Figures 2 (a) and (b) depict the 
heteronuclear PA spectra near the dissociation limits 
25i/2+4P3/2 and 2Si/2+4Pi/2 for PA detunings Apa be- 
tween and —325 GHz and between and —60 GHz, re- 
spectively. The graphs represent, respectively, an average 
of ~ 6 and ~ 20 recorded spectra for noise reduction and 
have been recorded in pieces and stitched together. The 
spectra contain 68 resonances whose contrasts decrease 
and whose mutual separations increase with increasing 
detuning. The maximum contrast amounts to ~ 35% 
and is obtained for a detuning of Apa = —14.4 GHz (see 
Fig. [2] (d)). The observed resonance widths (FWHM) 
vary between 80 and 300 MHz, primarily due to unre- 
solved molecular hyperfine structure. 

We have also recorded the heteronuclear -fhomonuclear 
PA spectra appearing on the K fluorescence signal, which 
contain all the resonances of fig. [2] (a,b) as well. A com- 
parison between the two spectra is shown for a small 
part in fig.[2](c). This figure shows comparable contrasts 
for the heteronuclear ^Li'*"K* and homonuclear PA 
signals. We identify the observed ^°K2 resonances as ex- 
citations to 0+ states [3l| . 

In the heteronuclear spectra of fig. [2] (a,b) we identify 
seven vibrational series (labeled with numbers), corre- 
sponding to seven of the eight molecular potentials disso- 
ciating to the 2S'i/2+4i'i/2,3/2 asymptotes (see fig.[l](b)). 
Each series contains up to five resonances, which appear 
in doublets or triplets due to resolved rotational struc- 
ture. This structure is shown more clearly in fig. [2] (c) 
for a particular vibrational state. Some of the observed 
rovibrational resonances have a further substructure re- 
sulting from hyperfine interactions, which is shown for a 
particular resonance in fig. [5] (d). 

In fig. [2] (a,b) and table H] we present an assignment 
of the observed resonances, which was obtained by the 
combination of different assignment rules: the first is the 
rotational progression law E,ot = By[J{J + l)-n'^, with 
J = + 1, ... for Hund's case c molecules [28| com- 
bined with our theoretical calculations of the rotational 



constants By. It allowed us to identify rotational progres- 
sions and to assign some J and based on the rotational 
spacing. The identification of the fl = 2 vibrational series 
(series 1 in fig. [5] (a)) is particularly easy, because only 
two rotational lines per vibrational level are expected, as 
opposed to three for all other series. The second is the 
semi-classical LeRoy-Bernstein (LRB) law |3^ Is^l (see 
Eq. ([1])) combined with the available calculated Cq co- 
efficients [2^ [l^]. It allowed us to identiiy vibrational 
progressions and to assign some v and based on the 
vibrational spacing. The third is the hyperfine structure 
law i?hfs oc n/[J(J -I- 1)] for = 1 and i^hfs ~ for 

= [33, HI] . It predicts small widths for resonances 
with ri = and particularly large widths for those with 
O = 1,J= 1 making their identification possible. The 
fourth is the expected similar contrast pattern of the ro- 
tational lines of the same vibrational series, which helped 
us to identify vibrational progressions. 

An application of the assignment rules allowed us to 
identify the observed vibrational series and to assign their 
quantum numbers except the parity a of the = elec- 
tronic states, a can be determined from an analysis of the 
relative strength of the rotational lines: due to the selec- 
tion rules, the parity of the total wave function of the sys- 
tem, i.e. the product of a and {—Y for the rotational part, 
changes sign during the transition. Further, a is con- 
served, namely only Ar^I]+(0+) 0+ and a^E+(0~) -> 
0~ are allowed for parallel transitions. In our experiment, 
s-wave collisions dominate, such that the total parity is 
-I- (— ) for the former (the latter) initial state. The par- 
allel transition X^T.+ {0+,£ = 0) ^ (0+,^ = 1) is thus 
allowed enhancing then the J = 1 line, while the paral- 
lel transition a^j:+{0-,e = 0) ^ (0-,€ = 1) is forbid- 
den. Under the same approximation, the perpendicular 
transition a^T,+ {l,i = 0) {0~,£ = 0) is allowed and 
enhances the J = line in the spectrum. Therefore we 
assign the $7 = series with pronounced (reduced) J = 1 
line to the excited 0+ (0~) states. 

IV. DISCUSSION 

Having assigned all observed resonances, the parame- 
ters of the different molecular potentials can be derived. 
We infer the Cg coefficients from the measured vibra- 
tional binding energies D — Ey = — (/lApA — i^rot) {D 
denoting the dissociation energy and Ey the energy of 
the vibrational level v), using the LRB formula (32l. Issj 

D-Ey^Ae{vD-vf, (1) 

with Aq = MTT^h^CQ/ [B{2/3,l/2)^/2iIC^f , where B 
denotes the Beta-function (5(2/3,1/2) w 2.587), ^ is 
the reduced mass, and vd the vibrational quantum num- 
ber at dissociation (a constant between and 1 so that 
the most weakly bound state has v = 1). Figure [3] shows 
the plots of the 1/3-rd power of the binding energies as 
a function of the vibrational quantum number for the 
five vibrational series dissociating to the 2S'i/2+4P3/2 
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-Ap^[GHz] 

FIG. 2: Heteronuclear PA trap loss spectra of ^Li^^K* below the asymptotes 25'i/2+4P3/2 (a) and 25'i/2+4Pi/2 (b). The 
spectra contain seven vibrational series (labeled = 1,...,7) with resolved rotational structure, whose assignment is given 
in table IT] (c) Zoom on the heteronuclear (upper trace, left axis) and heteronuclear+homonuclear (lower trace, right axis) 
PA spectrum below the 251/2+4^3/2 asymptote showing the rotational structure of the Q = — 3 vibrational state 

of ®Li*°K* (v denoting the vibrational quantum number counted from dissociation) and three vibrational Ou levels of '*°K2, 
which show a resolved hyperfine (* and o) but no rotational structure, (d) Zoom on the fl = = 2, J = 1 resonance of 

''Li'^'K*, showing a nearly resolved hyperfine structure. The PA detuning ApA is specified relative to the ^''K atomic transitions 
4Si/2(F = 9/2) ^ 4P3^2{F' = 11/2) (a,c,d) and 45i/2(F = 9/2) ^ 4P,^2{F' = 9/2) (b). 



asymptote. The plots are predicted to follow straight 
lines whose slopes yield: Cg = 9170 ± 940 a. u. and 
Ce = 9240 ± 960 a.u. for the dyad potentials n = 2, I^p, 
Cq = 25220 ± 600 a.u., Cg = 25454 ± 720 a.u. and 
Cq = 24310 ± 1710 a.u. for the upper triad potentials 
n = 1^^°™, 0+, 0- and Ce = 12860±660 a.u. for the lower 
triad potential 51 = 0+ (not shown in fig.El), respectively, 
where the uncertainties represent statistical uncertainties 
for the fits. These values are in good agreement with 
the respective theoretical values Cg — 9800 a.u., Cg — 
25500 a.u. and Cg = 13830 a.u. predicted by Bussery et 
al. [2^. The agreement with the values Cg — 9520 a.u., 
Cg = 22000 a.u. and Cg = 15420 a.u. predicted by Movre 
et al. [13] is not as good. The two predictions differ 
in their treatment of the interaction between the two 
asymptotes 2S'+4P and 2^+45", which is taken into ac- 
count in ref. (26j only, hinting its significance. 

The uncertainty of the derived Cg coefficients results 
from the following effects. First, the heteronuc lear na- 
ture of LiK and its small Cg coefficients lead to molecule 
formation at small internuclear separations (of i?cff — 
h/y/2jiB^ ^ 18 ao at Apa = -300 GHz) at which the 
exchange interaction and higher-order terms in the long- 
range multipole expansion of the molecular potential be- 
come important, which are neglected by the LRB law. 
Second, the small reduced mass of LiK leads to a low 
density of vibrational states and thus to a small number 
of states with long-range character available for fitting. 

The measured rotational splittings allow us to infer 
the rotational constants and to confirm the assignments 
above. They are shown in fig. U for the five vibrational 
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FIG. 3: Plot of the 1/3-rd power of the measured binding 
energies D — = — (/lApA — -Erot) (symbols) as a function 
of the vibrational quantum number counted from the disso- 
ciation limit for the five vibrational series dissociating to the 
25'i/2+4P3/2 asymptote. The slopes of the linear fits (solid 
lines) yield the dispersion coefficients Ce according to eq. ((l|. 
The nearly identical slopes of the triad and dyad potentials 
demonstrate the equality of the respective Cs coefficients. 

series below the 2S'i/2+4-P3/2 asymptote, together with 
their theoretical predictions, which we have derived from 
the potential curves of fig. [T] The agreement between 
the measured and predicted values is reasonable. The 
error bars account for the imprecision of the wavelength 
determination and of the resonance positions due to the 
unresolved hyperfine structure. Deviations from the the- 
oretical predictions are likely to be due to the multichan- 
nel character of the vibrational levels. 
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TABLE I: PA resonances of ''Li^^K* observed below the 
2<S'i/2±4Pi/2,3/2 asymptotes and their contrasts. A'^ denotes 
the number of the vibrational series given in fig. [2](a,b). 

We have determined the ^Li^^K* molecule formation 
rate from the steady-state depletion of the ®Li atom num- 
ber induced by PA. For the resonance shown in fig. [2] (d) 
we obtain a lower bound of /Spa^k-^l/^ ~ 3.5 x 10'' s^^ 
and a PA rate coefficient /3pa = (2.2± 1.1) x lO^^^ cmS/s. 
This coefficient is larger by about a factor of two than 
the one found in the experiment with RbCs* [l^, show- 
ing that PA rates for LiK* are much more favorable than 
previously expected (2l| . confirming the trend discussed 
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FIG. 4: Measured rotational constants (symbols) for the 
observed excited molecular states below the 2S'i/2±4P3/2 
asymptote and their theoretical predictions for computed vi- 
brational levels (dots, the lines serve to guide the eye), derived 
from the potential curves of fig. [T] 

in ref. [3S]. Using the approach described in ref. [s^] we 
estimate the total PA rate coefficient for a computed 1"p 
level with -21 GHz detuning to 1.6 x IQ-^^cm^/s [IH, 
in agreement with our measured value. The associated 
^Li'^^K* molecule formation rate is also found to be com- 
parable to that for "^^Kj, which is derived from fig. [2](c) 
to be ^ 5.3 X 10^ s^^, despite the much longer range of 
the excited "'"Kj molecular potential. Using our model, 
which reproduces the observed rovibrational structure, 
we infer the rates for decay into bound levels of the 
Ar^E+ state. The rates are found largest for the decay 
into the most weakly bound level, but are still signifi- 
cant for the decay into deeply bound levels such as the 
X^'S'^{v' = 3) level {v' counted from the potential bot- 
tom) for which it is 5 x 10"*^ s~^ (see fig. [T]). Since the 
^Li'^'^K* molecule formation rate saturates in our experi- 
ment at moderate PA intensities of ~ 80 W/cm^, efficient 
coherent multi-photon population transfers to the molec- 
ular rovibrational ground state can be expected. 

V. CONCLUSION 

In summary, we have investigated single-photon pho- 
toassociation of excited heteronuclear ^Li^'^K* molecules. 
We have recorded photoassociation spectra and assigned 
all observed resonances. We have derived the long- 
range dispersion coefficients and rotational constants, 
which agree with the theoretical predictions of ref. [2^ 
and our calculations, respectively. In particular, we 
have observed large formation rates for the heteronuclear 
^Li'*°K* molecules which are comparable to those found 
for homonuclear ^^Kj. These rates promise efficient cre- 
ation of rovibrational ground-state molecules and show 
that photoassociation is an attractive alternative to Fes- 
hbach resonances, since those have a ver y sm all width for 
^Li^°K and are thus difficult to control [ij]. 

For future research it will be interesting to combine 
ours and previously recorded data on LiK and to refine 
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the molecular potentials. Sp ectroscopic data is avail- 
able for the potential l^II [Sa], which correlates with the 
— 1"P potential, for which we measured the binding 
energies of the five previously undetermined least-bound 
vibrational states, such that a complete set of vibrational 
levels is now available for a high-precision refinement of 
this potential. 
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